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We report smooth relativistic Hartree-Fock pseudopotentials (also known as averaged relativistic 
effective potentials or AREPs) and spin-orbit operators for the atoms H to Ba and Lu to Hg. 
We remove the unphysical extremely non-local behaviour resulting from the exchange interaction 
in a controlled manner, and represent the resulting pseudopotentials in an analytic form suitable 
for use within standard quantum chemistry codes. These pseudopotentials are suitable for use 
within Hartree-Fock and correlated wave function methods, including diffusion quantum Monte 
Carlo calculations. 



PACS numbers: 71.15.Dx, 31.15.-p, 02.70.Ss 

Pseudopotentials or effective core potentials (ECPs) 
are commonly used within electronic structure calcula- 
tions to replace the chemically inert core electrons. The 
influence of the core on the valence electrons is then 
described by an angular-momentum-dependent effective 
potential, leading to greatly improved computational effi- 
ciency in ab initio calculations for heavy atoms. The use 
of pseudopotentials is well established within Hartree- 
Fock (HF) and Density Functional Theory (DFT), and 
in correlated wave function calculations. 

Our main interest is in diffusion quantum Monte Carlo 
(DMC) calculations pi* 2 - This technique provides an ac- 
curate solution of the interacting electron problem for 
which the computational effort scales with the number 
of electrons, N, as approximately N 3 , which is better 
than other correlated wave function approaches. Un- 
fortunately, scaling with atomic number, Z, is approx- 
imately^ Z 5 ~ 6 - 5 . The use of a pseudopotential reduces 
the effective value of Z, making DMC calculations feasi- 
ble for heavy atoms. 

There is evidence that HF pseudopotentials give bet- 
ter results within DMC than DFT pseudopotentials; 5 It 
appears that the complete neglect of core-valence corre- 
lation within HF theory leads to better pseudopotentials 
than the description of core-valence correlation provided 
by DFT. Moreover, core-valence correlation can be in- 
cluded within correlated wave function calculations per- 
formed with HF pseudopotentials by using core polariza- 
tion potentials£i£i£ Core polarization potentials mimic 
the effects of dynamical polarization of the core by the 
valence electrons, as well as static polarization effects 
due to the other ions. We would therefore like to use 
HF pseudopotentials in our DMC calculations, prefer- 
ably constructed from Dirac-Fock (DF) theory in order 
to include the relativistic effects which are significant for 
heavy atoms. 

Standard quantum chemistry packages are convenient 
for generating the "guiding wave functions" required in 
DMC calculations. We would therefore like our pseu- 
dopotentials to be available in the standard parameter- 
ized form of a sum of Gaussian functions multiplied by 
powers of the electron-nucleus separation. 

Extensive sets of parameterized pseudopotentials are 



available in the literature, but they have generally been 
constructed with different goals to ours. Relativistic 
pseudopotentialsS. generated within DFT and the lo- 
cal density approximation have been available for some 
timei^ but, as mentioned above, it appears that HF pseu- 
dopotentials are superior for our purposes. 

Pseudopotentials defined within HF theory have been 
published by Hay and Wadtii for much of the periodic 
table, and Christiansen and coworker s 12 : 13 have gener- 
ated similar HF pseudopotentials from DF atomic cal- 
culations, thereby including relativistic effects. "Energy 
consistent" HF pseudopotentials including relativistic ef- 
fects have also been developed by the Stuttgart and Bonn 
groups^, and are also publicly available^ There is con- 
siderable freedom in constructing pseudopotentials, and 
many of those available in the quantum chemistry liter- 
ature diverge at the nucleus, normally as 1/r 2 or 1/r. 
This singular behaviour leads to large "time-step" errors 
and even instabilities in DMC calculations^ Moreover, 
we cannot imagine that this singular behaviour is advan- 
tageous in quantum chemistry methods, as it leads to 
pseudo wave functions with behaviour at small r which 
cannot be described within a Gaussian basis set. All 
of the HF pseudopotentials discussed above possess such 
singularities at the nuclear site. 

Greeff and Lester— and Ovcharenko et aZ,~ have gen- 
erated non-singular HF pseudopotentials for the atoms 
B-Ne and Al-Ar, which are available in a suitable pa- 
rameterized form. These are constructed explicitly for 
use in DMC calculations, and satisfy some of our crite- 
ria. Unfortunately these pseudopotentials cover only a 
small part of the periodic table, do not have d angular 
momentum channels, do not include relativistic effects, 
and are not particularly smooth even though they are 
non-singular at the origin. 

In this paper we report the generation of a library of 
accurate HF pseudopotentials which are non-singular at 
the origin, include relativistic effects, are parameterized 
in a form appropriate for quantum chemistry packages, 
and are smooth to aid transferability. We also require our 
pseudopotentials to be well-localized (by which we mean 
non-local only in a small region around the nucleus), as 
the evaluation of the non-local energy within DMC is 
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expensive. 

We also seek to quantify the errors present in our pseu- 
dopotentials. Transferability errors are expected to be 
the most significant, but errors due to other sources are 
also assessed. These errors result from the approxima- 
tions required to construct HF pseudopotentials from all- 
electron (AE) DF data, from the removal of the unphys- 
ical long ranged non-local behaviour discussed by Trail 
and Needsil, and from the parameterisation of the pseu- 
dopotentials. 

The rest of this paper is organized as follows. In sec- 
tion [T] we describe the form of the relativistic pseudopo- 
tentials obtained from AE DF orbitals by inversion of 
the HF equations. In section [II] we describe the partic- 
ular method we use to generate pseudopotentials. The 
unphysical non-local tails of these pseudopotentials are 
removed and we parameterize them accurately. In sec- 
tion IIIII we present and analyse pseudopotentials gener- 
ated for the atoms H-Ba and Lu-Hg, and present the 
results of test calculations for a selection of these atoms, 
making comparisons with AE DF results and with re- 
sults obtained using the pseudopotentials of Christiansen 
et al. We also generate parameterized spin-orbit (SO) 
pseudopotentials and compare results of calculations us- 
ing these with AE DF results, and with results obtained 
using the SO pseudopotentials of Christiansen et al. We 
draw our conclusions in section HVl 

Atomic units are used throughout, unless otherwise in- 
dicated. 



FORM OF THE PSEUDOPOTENTIALS 



The coupled radial DF equation o 18 i 19 are 
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where a — 1/c is the fine structure constant, i is the 
state index, and Gi and Fi are the radial components of 
the major and minor (or large and small) parts of the 
Dirac orbitals. The i index is unique to each orbital, and 
hence the non-zero integer quantum number ki can be 
associated with each i, 



h = k + \ ki>0 
-(*< + !) = -& + §) h<0. 



(3) 



The orbital dependent effective potentials, V? and 
V^f' F , are functional of the set {G,F}. 

In what follows we also consider the standard approx- 
imation to Eqs. (11I2[) for valence orbit alsi^^ Taking q 
and Vi F to be small in Eqs. (|1I2|) . as is the case for 
the valence orbitals outside of the core region (this can 
be taken as one of the criteria for the core/valence par- 
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Equations (|4]) and ([5]) are equivalent to Eqs. (|H2p up to, 
but not including, terms of order a 2 , and it should be 
noted that a self-consistent eigenstate of Eqs. (|H2p does 
not satisfy Eqs. (|4l5p exactly in any region of space. The 
effective potential in Eqs. @ and Q is given by 
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where the first term arises from the nuclear charge and 
the second is the Hartree potential due to the total elec- 
tron density, p. The third term is the effective exchange 
potential acting on the major part of the radial Dirac 
orbital. Both the Hartree and exchange terms are func- 
tional of {G,F}, and the effective exchange potential is 
different for each orbital (the cancelling self-interaction 
has been included in the exchange and Hartree terms). 

To generate a non-relativistic HF pseudopotential from 
the AE DF solutions of Eqs. (|H2p we follow the proce- 
dure described by Kleinman^ The AE DF atomic solu- 
tions are partitioned into core states whose influence is 
represented by the pseudopotential, and valence states 
which are represented by the pseudo-orbitals. We then 
construct non-relativistic pseudo-orbitals which preserve 
desired properties of the original AE DF valence orbitals. 

The pseudo-orbitals are chosen to satisfy the 
Schrodinger-like equation 
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where the pseudo-orbital, ^>j, is a scalar orbital and the 
corresponding eigenvalue ej is equal to that of Eqs. (|H2p 
for the AE orbital i. We may relate Eq. to the Dirac- 
Fock equation by taking the limit a — > in Eqs. (|4l5p . 
which leads to a scaler Schrodinger equation. 

The effective potential, for state i is then given 

by 
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where pP seudo is the charge density obtained from the 
occupied pseudo states. This defines the pseudopoten- 
tial for orbital i, Vi(r), in terms of the pseudo-orbital 
</>,:. Given <pi and ej, Vi(r) may be obtained by direct 
inversion of Eqs. (O and ([8|). 

In the non-relativistic case we would construct valence 
pseudo-orbitals, </>,;, which are equal to the equivalent AE 
orbitals outside of a 'core radius', r ci . In the relativistic 
case this procedure must be modified, since the DF or- 
bitals consist of two components. We could simply take 
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(f>i = Gi for r > r C i, but this is unsatisfactory because it 
implies that charge has been removed from the system. 
Instead we take the form 
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[Gi(r) 2 + Fi(r) 2 ]* r > r c 
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which preserves the valence charge density outside of the 
core region. Inside r c j the orbitals are given by the, as 
yet, unspecified function /j. This function is generally 
chosen such that <pi is node-less, smooth at r C i to a certain 
order of differentiation, and that it satisfies the norm- 
conservation condition, that is, the total charge inside of 
r C i is the same as for the AE orbital. 

The pseudopotentials Vi are not quite appropriate for 
use in non-relativistic calculations as they provide dif- 
ferent pseudopotentials for the same I values. However, 
the valence states chosen to construct the pseudopoten- 
tial have unique k numbers, and therefore the index k 
is interchangeable with the pseudo-atom orbital index, 
i. Taking the average of the pseudopotentials with the 
same I quantum numbers weighted by the different j de- 
generacies, we obtain the averaged relativistic effective 
potential (AREP^, 
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This includes all of the Dirac relativistic effects except the 
SO coupling. The SO operator, Vf , may be expressed 
in terms of a SO pseudopotential, V, s °, and angular mo- 
mentum operators, 
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To use the pseudopotential in a calculation for a 
molecule or solid it is expressed in terms of projection 
operators, and separated into local and non-local parts, 

Vpseudo — ^oca;( r ) 

+ EE \ Y i™){Vnr)-V? ocal {r))(Yim\. 
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Orbitals with / > l max feel the local potential, Vf ocal . 



II. GENERATION OF THE 
PSEUDOPOTENTIALS 

We take the core states to be the 'noble core' (for ex- 
ample, for Si we take Is2s2p as core states and 3s3p as 
valence states) or noble core plus a filled d shell. We gen- 
erate pseudopotentials for the s,p and d channels. We use 
the atomic ground state configurations to obtain chan- 
nels for which the corresponding AE valence orbitals are 



bound. For channels which have no corresponding AE 
valence orbitals in the ground state we use the excited 
state configurations of Bachelet et al*£, which provide 
an appropriate bound state. Fischer's interpretation of 
fractional occupation numbers within HF and DF the- 
ory^ is used. In most cases we choose the core radii, r c i, 
to be 0.9 of the distance from the outermost node to the 
outermost maximum, although for some cases we choose 
smaller radii. 

To define our pseudo-orbitals we use the Troullicr- 
Martins scheme 22 where the pseudo-orbitals within the 
core region are given by 
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and the coefficients C2 m are to be determined. The fac- 
tor r li+1 ensures that no 1/r 2 singularity occurs in the 
effective potential (and resulting pseudopotential), and 
Ci = ensures that no 1/r singularity is present. All of 
the C2m+i terms have been excluded to prevent the ap- 
pearance of a cusp of any order in the pseudopotentials 
at the origin, so improving the asymptotic behaviour of 
the pseudo-orbitals and pseudopotentials in momentum 
space. The seven coefficients in Eq. (Tl4"|) are determined 
by the conditions of: 



1 . norm-conservation within r r 
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2. continuity of 4>i(r) and its first four derivatives at 

' ci j 

3. zero curvature of the screened potential at the ori- 
gin, 
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In the core region we obtain by inversion of Eq. J7]| , 
and obtain V% by 'unscreening' using the pseudo- 

orbitals. 

Outside of the core region we could in principle ap- 
ply the same inversion procedure, but this would involve 
taking numerical differences of small quantities, which is 
prone to errors. Instead we note that the eigenvalues of 
the AE DF and pseudo valence states are equal, and that 
we expect Eqs. (I4I5P to be an extremely good approxima- 
tion for valence states in the region outside of the core. 
We may then use Eqs. (g]), ©, Q, and © to obtain 



Vi = +V h [p-p 



pseudoi 



V^[{G,F},ki}Gi v x G l{<f>,Q},h)<t>i 



G, 



~G 



(17) 



4 



for r > r c i. The final term in this expression arises from 
the kinetic energy, and we neglect this, which introduces 
an error of order a 2 . This approximation may seem un- 
necessary, but errors of this order are already present 
if we assume that the pseudo-valence electrons may be 
treated as scalar relativistic, and it has the advantage 
that the valence charge density outside of the core region 
is preserved. We also found that this error, when de- 
tectable, is smaller than that due to the fact that Eq. ((4]) 
is already an approximation, even outside of the core re- 
gion. Some authors^ have chosen fa = XGi outside of 
the core, so that the kinetic terms in Eq. (jTTJ) are zero. 
The A constant may be chosen to ensure that fa is nor- 
malized. This would cause a negligible (in most cases 
undetectable) change in the pseudopotentials generated 
here. 

All DF and HF calculations were carried out using fine 
radial grids and therefore do not suffer from basis set er- 
rors. We used the DF code of Ankudinov et al££ for the 
AE DF calculations, the HF code of Fischer— for the AE 
HF calculations, and the pseudopotential HF calculations 
were performed using our own code. Grid and conver- 
gence parameters were chosen so as to achieve 10-digit 
accuracy in the total energy, and the ground state ener- 
gies were checked against those of Visscher and Dyall 2 ^, 
and Librelon and Jorge^ The above procedure provides 
us with relativistic screened pseudopotentials for each of 
the valence states considered and, from Eq. (fT7|) . the V[. 
We then use Eq. flTUJ) and |T5J to obtain the AREP and 
SO potentials. 



A. Non-local asymptotic behaviour 

In a previous paper— we examined pseudopotentials 
constructed within HF theory using this inversion pro- 
cedure, and found them to be non-local over all space, 
the deviation from the ionic Coulomb potential remain- 
ing finite as r — > oo. As the presence of this 'extreme 
non-locality' is a consequence of the non-locality of the 
exchange interaction in HF theory, and the DF exchange 
interaction is of essentially the same form, this effect also 
occurs for the pseudopotentials defined here. 

The presence of this long-ranged non-locality leads to 
a loss of transferability and problems in defining a total 
energy for extended systems. Here we remove the 'ex- 
treme non-locality' using the method presented in our 
previous paper—, which we briefly summarize here. 

In order to generate a new version of the pseudopoten- 
tial, Vj P ' (oc (r), which is non-local only close to the atomic 
core, the original pseudopotential is transformed using 

+Vf(r) r < r c 

(li(r) + V l p (r)~V h [p core } + ^) 

+ Vh[pcore] - f r > r c , 

(18) 



where Vh[p CO re\ — Z/r is the ionic potential and r\~^ — 
max; [r c i\ /16 is a parameter which specifies the length 
scale over which the transformed pseudopotential be- 
comes local outside of the core region. The function 7; is 
given by 

, , / Qi + Pir 4 (l - 3^r 2 ) r<r c 
li{r) = < r4 V sr ° J (19) 

I Ql + Pi -f r>r c , 

where the parameters qi and pi are chosen such that the 
HF eigenvalues of the original pseudopotential are pre- 
served, and that the logarithmic derivatives of the orig- 
inal pseudo-orbitals at r C i are preserved to high accu- 
racy?^ Although the norm of the original pseudopoten- 
tial is not exactly preserved this was not found to signif- 
icantly affect transferability. 

B. Expansion in a Gaussian basis 

For applications within quantum chemistry codes it is 
normally necessary to have the pseudopotentials avail- 
able as an expansion in Gaussian functions. The pseu- 
dopotentials described so far are tabulated on a radial 
grid, and in this section we develop a Gaussian fitting 
procedure which is both accurate and reliable. 

Our aim is to find an expansion of the localized AREP 
yp,ioc j n a g ^ anc [ arc j Gaussian form, 

yP _ V" A , r n ql -a ql r 2 _ / Z val /r + I = focal 

q ^ £ Local ' 

where Z va i is the number of valence electrons in the neu- 
tral atom. A set of n q \ is chosen such that the pseudopo- 
tential can be represented accurately by the expansion, 
with no singularity present at the origin. Once we have 
this expansion it may be applied in the same projector 
form as Eq. Q13p . It is important to note that we ex- 
pand the localized tabulated pseudopotential in terms of 
Gaussian functions. This is desirable since the removal 
of the long-ranged non-local tail significantly improves 
the accuracy and stability of any fitting procedure using 
localized basis functions, such as the Gaussian expansion. 

In order arrive at an algorithm to generate our expan- 
sion we must decide what properties a 'good' Gaussian 
expansion must have. We require a good expansion to 
reproduce both the eigenvalues and orbitals of the tabu- 
lated pseudopotential to high accuracy, and so reproduce 
the scattering properties of the tabulated pseudopoten- 
tial. We require a good expansion to be close to the tabu- 
lated pseudopotentials in the least-squares sense, because 
such deviations tend to make the fitted ones less smooth 
and reduce the transferability. We also require the pa- 
rameterized pseudopotentials to be non-local only in a 
region near the core, although this criterion is partly in- 
cluded in the previous two, together with the use of the 
localized tabulated pseudopotential. These are stringent 
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requirements and we therefore require a rather larger ex- 
pansion than has been used in earlier work. Finally we 
require an expansion which is compatible with quantum 
chemistry codes, and we have successfully tested our pa- 
rameterized form in both the CRYSTAL 2 ! and GAUSSIAN^ 
codes. 

We chose the local part of the pseudopotential to be 
I = 2 for all atoms, and to obtain an accurate (in the least 
squares sense) representation of the pseudopotential we 
chose qmax = 8, with n q i = —1,0, 1,2 for the local part 
and n q i = 0, 1,2 for the non-local parts. This gives 16 
parameters for each channel. It is helpful to reduce this 
number by imposing constraints on the functional form of 
the pseudopotential. We require that is finite and has 
zero derivative at the nucleus, as did Greeff and Lester^ 
In addition we require the second derivative of to be 
zero at the nucleus and that VJ p (0) = Vf' oc (0). The lat- 
ter two conditions make smoother because they force 
it to be closer to the original tabulated pseudopotential. 
These conditions are imposed by reducing the number of 
free A q i parameters. 

Since Eq. ([20]) is linear in the {A q {\ it is also possible 
to use Singular Value Decomposition (SVD) to define the 
remaining A q i parameters as those which give the mini- 
mum least squares deviation from the original tabulated 
pseudopotential. 

A further (nonlinear) least squares minimization with 
respect to the remaining free parameters {a q i} provides a 
'good fit' to the tabulated pseudopotential, but does not 
provide an expansion which accurately reproduces the 
eigenstates of the original pseudopotential. We take the 
{a q i} parameters obtained from this least squares proce- 
dure as a starting point for a second stage of optimiza- 
tion. Here we employ a generalization of an algorithm 
developed by Barthelat et al~ 

To progress further we perform a HF atomic calcula- 
tion with LS coupling and using the tabulated pseudopo- 
tential Vf' oc , yielding the pseudo-states {<$ ,ef }. For a 
given set of parameters for each channel, we also perform 
a similar HF atomic calculation using the parameterized 
pseudopotential Vf to give the pseudo-states {0f ,ef }. 
We then define a functional of the pseudo-states of both 
the tabulated and parameterized pseudopotentials, 



where 



Oi 



efl#X#l-ef|tfW 



(21) 



(22) 



which is a function of the {a q i} for all I which corre- 
spond to valence states in the atomic configuration con- 
sidered. We vary the available a q i of the parameterized 
pseudopotential to search for a minimum of £ where the 
overlap (0f |</>f ) is maximized and the error in the eigen- 
value |ef — ef | is minimized separately for each available 
I. In the original formulation of Barthelat et al^ no 



sum over states was present in Eq. (|21l) . and each chan- 
nel was optimized separately. We found that optimizing 
each channel separately did not allow convergence to be 
achieved to a high enough accuracy for the majority of 
atoms, whereas optimization over all states provided re- 
liable convergence. 

A standard quasi-Newton minimization algorithm was 
used in the optimization, and we found that the stability 
and efficiency could be improved by introducing 
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and minimizing with respect to the {a! ql \. 

We considered the optimization to be successful when 
the pseudo-states of the parameterized pseudopotentials 
satisfied the conditions 
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for the transition metal atoms. 

On implementation it became apparent that, for most 
atoms, the final values of the parameters are surprisingly 
sensitive to their initial values, suggesting that £ has 
many minima. We also found that, without the SVD def- 
inition of the A q i parameters, the rescaling of Eq. (j2"3"|) 
or the concomitant optimization implied by the sum in 
Eq. (|21[) . this procedure did not converge for many atoms, 
suggesting that these aspects of the algorithm are nec- 
essary to provide a functional with minima which are 
distinct and well defined. The presence of many locally 
optimum parameter sets is also an advantage, as when 
a set of optimized parameters does not satisfy our crite- 
ria a new contender may easily be obtained by making a 
small change to the initial parameter set. 

As an example we consider C, for which our param- 
eterized pseudopotential is shown in Fig. [TJ The s and 
p channels were obtained from the neutral ground state, 
while the d channel was obtained from an excited ionic 
state (ls 2 2s a75 2p 1 00 3d - 25 ) taken from Bachelet et al. 10 
Inset in the same figure is the difference between the 
parameterized pseudopotential and the original TM (ex- 
tremely non-local) pseudopotential. This figure demon- 
strates that our parameterisation accurately reproduces 
the original pseudopotential, with the greatest difference 
occurring at small r. Most of this difference is not due 
to the fitting procedure - it is due to the change in the 
pseudopotential required to enforce locality outside of a 
small radius (in this case r; oc = 1.2 a.u.) while preserving 
the eigenvalues and logarithmic derivatives of the original 
pseudo-states. 

In Fig.[2]the s and p components of our C pseudopoten- 
tial are compared with the non-singular (non-relativistic) 
pseudopotential published by Ovcharenko et al} 6 Fig. [3] 
shows that our C pseudopotential is smoother than that 
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r (a.u.) 



FIG. 1: Our parameterized pseudopotential for C. On the 
scale shown it is indistinguishable from both the original tabu- 
lated TM pseudopotential and its localized version. The inset 
shows the difference between the parameterized and original 
TM pseudopotentials. 



of Ovcharenko et al, and we found this to be the case 
for all of their pseudopotentials. This is probably due to 
the extra constraints we apply in our construction of the 
pseudo-orbitals . 

Fig. [3] shows the s pseudo-orbitals resulting from our 
C pseudopotential, together with those from the pseu- 
dopotentials of Ovcharenko et al*ZL (OAL), Christiansen 
et alH (PC), and those of the Stuttgart group^ . (ISP). 
The +l/r 2 divergence of the PC pseudopotential forces 
the s pseudo-orbital to go rapidly to zero at the nucleus, 
while the — 1/r divergence of the ISP pseudopotential 
forces the s pseudo-orbital to have a cusp at the nucleus. 
Our pseudopotential gives the smoothest pseudo-orbital 
of those plotted. 

We have generated pseudopotentials using the scheme 
described above for the atoms H to Ba and Lu to Hg. We 
include H and He even though they have no core electrons 
because the smoothness of the resulting pseudopotentials 
may make them useful in some circumstances. 

To provide a quantitative measure of the errors intro- 
duced by the various approximations used in the con- 
struction of the pseudopotentials we investigated how 
well the pseudopotentials reproduce the AE atomic re- 
sults at the various levels of approximation. Atomic LS 
and jj coupled HF calculations were performed with the 
pseudopotentials using our own code. Here we report de- 
tailed results for the Li, Be, C, Si, Ti, Cr, Fe, Br, Mo, 
Ag, and Sb atoms, which are representative of the gen- 
eral accuracy achieved. Note that we have deliberately 
chosen several transition metal atoms as they generally 
exhibit the largest errors. 

If there were no approximations in our pseudopoten- 
tial construction procedure, then a HF pseudopotential 



FIG. 2: Comparison of the s and p channels of the parame- 
terized C pseudopotentials generated within this paper (Our) 
and by Ovcharenko et aZ.— (OAL). (The OAL pseudopoten- 
tial does not have a d channel for us to compare with.) 



calculation with the exchange interaction defined using 
jj coupling should exactly reproduce the AE DF eigen- 
values of the occupied levels of the ground state. In 
practice, however, this is not achieved, and the pseu- 
dopotential generation method we use may be viewed as 
introducing four types of error. These are shown in Ta- 
ble U together with the AE DF eigenvalues, and the jj 
coupled HF eigenvalues resulting from the final param- 
eterized pseudopotential. For both the errors and the 
eigenvalues themselves we take an average over values 
corresponding to the same I weighted by the different 
j-degeneracies, just as we did to construct the AREP. 

The first type of error, Aei, is introduced by the fact 
that the pseudo-orbital is not an exact solution of Eq. ^ 
outside of the core region. The second error, Ae2, is in- 
troduced by the j-averaging of the pseudopotential. This 
error results from the exclusion of SO coupling from the 
pseudopotential, and is due to the higher order effects 
of SO coupling. The third error Ae3 is introduced by 
the localization of the pseudopotential (by construction 
this error would be zero for LS coupling). Finally, the 
fourth error, (Ae^), is introduced by the imperfect pa- 
rameterisation. These four errors are shown in Table U 
and are defined such that their sum is the total error 
of the (averaged) eigenvalues resulting from the parame- 
terized pseudopotential. Results are given for the atoms 
listed above, together with the separate averages over the 
s, p and d block atoms considered in this paper. 

The total error is small for all atoms, it is negligible for 
the s-block atoms and is largest for the heaviest transi- 
tion metal atoms. Generally, the largest total errors are 
dominated by the error due to the removal of SO cou- 
pling, this error increasing rapidly with increasing atomic 
number (compare Ae2 for Mo and Ag). Aei + Ae3 + Ae± 
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FIG. 3: Comparison of the s pseudo-orbital from the neutral 
ground state of our parameterized C pseudopotential with 
those generated from the parameterized pseudopotentials of 
Ovcharenko et al^A (OAL), Pacios and Christiansen^^ (PC) 
and Igel-Mann et aZ.— 

remains below 1CP 4 a.u., with the exception of a few 
heavy transition metal atoms such as Mo. Errors due to 
the approximation of the Dirac equation as a Schrodinger 
equation in the valence region remain consistently small, 
and errors due to localization and parameterisation are 
largest for the transition metal atoms. 

In Table |TT] we compare the j-averaged eigenvalues re- 
sulting from AE DF calculations, eigenvalues resulting 
from HF calculations using our soft parameterized pseu- 
dopotentials, and eigenvalues resulting from the singular 
PC pseudopotentials. For all comparisons with PC pseu- 
dopotentials we choose those constructed using the same 
valence/core partition as our own. The j-averaged va- 
lence eigenvalues of the original AE DF atom are well 
reproduced by both the soft and singular pseudopoten- 
tials, with our soft pseudopotentials performing slightly 
better. We conclude that our pseudopotentials have been 
parameterized successfully. 



III. TESTING OF THE PSEUDOPOTENTIALS 

The analysis in the previous section demonstrates 
that our pseudopotentials have been parameterized accu- 
rately. Of course this does not demonstrate their trans- 
ferability, and it is this issue which is addressed in this 
section. 

Although designed to reproduce jj coupling results, 
our pseudopotentials are intended for use within LS 
coupling calculations, and therefore we have calculated 
HF excitation energies and ionization energies within 
LS coupling for both our pseudopotentials and those of 
Christiansen et alr^ (PC). Both of these types of pseu- 



dopotential contain relativistic effects at the DF level, 
and therefore we compare these results with AE HF re- 
sults in which scalar relativistic effects are incorporated 
perturbatively. The calculations were performed using 
Fischer's code with the Brcit-Pauli scalar relativistic cor- 
rection added to the total energy, correct to order a 2 . 
This order of correction would not be particularly accu- 
rate for AE calculations on the heavier atoms, but these 
errors arise largely from the deep core electrons and do 
not contribute significantly to the energy differences be- 
tween the low- lying excited states considered here. 

In Tabic IIIII we compare excitation energies for the 
same atoms considered in Tables U and [TTJ Excitation 
energies obtained for the AE atoms are compared with 
those from our soft AREP parameterized pseudopoten- 
tials and the singular PC pseudopotentials. We find that 
our pseudopotentials reproduce the excitation energies to 
the same accuracy as those of PC. 

We also consider the errors for our full set of pseu- 
dopotentials in reproducing the AE DF excitation ener- 
gies. It is apparent that the excitation energies are repro- 
duced most accurately for the s-block atoms with an av- 
erage absolute error of 0.00084 a.u. (0.022 eV), that they 
are less well reproduced for p-block atoms with an aver- 
age absolute error of 0.00344 a.u. (0.090 eV), and that 
the largest errors occur for the transition metal atoms, 
with an average absolute error of 0.01204 a.u. (0.32 eV). 
The large error value for the transition metals is mainly 
due to the 5d transitions - the average absolute errors 
are 0.00687 a.u. (0.18 eV), 0.00754 a.u. (0.20 eV), and 
0.02116 a.u. (0.55 eV) for the 3d, id, and 5d transitions 
metals respectively. Further examination of the results 
reveals that these errors are inherent in the pseudopo- 
tentials themselves, and not due to the errors examined 
in the previous section. Our pseudopotentials and the 
PC pseudopotentials also appear to produce very simi- 
lar errors, providing support for this interpretation. For 
the transition metal atoms it should be possible to re- 
duce these errors by including the 3 s and 3p electrons as 
valence states, giving "small core" pseudopotentials.-2i 

In Table HVl we compare the first ionization potentials 
calculated for the atoms considered in the previous ta- 
bles. The results are of similar quality to those for the 
excitation energies given in Table HTT1 and show a similar 
trend that the errors tend to increase with Z. The PC 
pseudopotentials perform slightly better for the ioniza- 
tion potentials than our soft pseudopotentials, but the 
difference does not appear to be significant considering 
the number of atoms studied. 

Our pseudopotential generation procedure also pro- 
vides a description of SO splitting in terms of the SO 
pseudopotential, V^ so defined in Eq. (Tj"2"]) . We may define 
a Gaussian parameterisation of this as 

Vf = B ql r n -»e- b -» r2 , (26) 
i 

where Vf° is the fitted potential and the B q i and b q i 
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play the same role as the A q i and a q i in the AREP. No 
Coulomb potential terms are required, and no local po- 
tential term is required since Vf° — for I > l max is 
a physically reasonable approximation. In addition we 
require B q i — for I — 0, since — for I = 0. 

It would be possible to follow the same route of local- 
ization and parameterisation used for the AREP. This 
would require jj coupled HF calculations to be used in 
the localization and fitting procedures of section[H](as op- 
posed to LS coupled HF), so we could seek a which 
reproduces the SO splitting of the original AE DF cal- 
culation. Here we use a simpler procedure and apply 
least squares fitting with respect to {B q i,b q i}, with the 
initial {b q i} taken as the values of the equivalent AREP 
parameters. 

We considered the selection of atoms and ions shown 
in Table |V] The SO splitting is defined as the differ- 
ence between the energies of states where J differs by 
one due to the transfer of an electron between orbitals of 
the same I, but different j quantum number. The differ- 
ence in energy of the states due to the SO interaction was 
obtained from AE DF calculations (not including further 
relativistic corrections), and from jj-coupled HF calcu- 
lations with the parameterized pseudopotentials. For the 
pseudopotential calculations V p and Vf° are used to con- 
struct a Relativistic Effective Potential (REP) that de- 
pends on j - the parameterized equivalent of Vi in section 
U - and this REP provides the external potential. Since 
the SO splitting is not obtained perturbatively any dif- 
ferences between AE and AREP+SO results are solely 
due to the REP. 

We calculated energy differences between 2P 1 / 2 and 
2P 3 / 2 states for s and p-block atoms, and between 2D 3 / 2 
and 2D 5 /2 states for ci-block atoms. In Table Ivl we com- 
pare the SO splitting energies obtained from AE DF 
calculations, those obtained using our parameterisation 
of V p and Vf , and those obtained using the (singular) 
parameterized SO pseudopotentials published by Pacios 
and Christianseni 13 i 32 ' 33 Our AE results differ somewhat 
from those reported by Christiansen et al£^, who used 
the code of Desclaux 3 -. These discrepancies appear to be 
due to the finer radial grids and higher tolerances used 
in our calculations. We have also found discrepancies be- 
tween our results and Christiansen's et aZ. 34 using their 
REPs, which we believe is due to the same effect. We 
also note that our code reproduces the accurate ground 
state energies of Visscher and Dyall^, but the code of 
Desclaux 3 ^ does not. 

We find that in most cases our pseudopotentials repro- 
duce the DF splitting energies more accurately than those 
of PC, with an average absolute percentage error of 4.8 % 
compared to 8.2 % for PC. The largest percentage error 



is for Li (27 %, compared with 34 % from PC), where the 
splitting is small. For the rest of the atoms considered, 
the errors are < 11 %, compared with < 24 % for PC. 

Parameter values for our pseudopotentials and 
SO pseudopotentials are provided as supplementary 
material 3 ^, and are given to the same numerical preci- 
sion used to generate the results in this paper. 



IV. CONCLUSIONS 

We have developed smooth HF pseudopotentials which 
emulate the influence of relativistic core electrons for the 
atoms H to Ba and Lu to Hg. These are the Averaged 
Relativistic Effective Potentials and SO pseudopoten- 
tials of Kleinman^ and Bachelet and Schliiterj2£ We use 
the Troullier-Martins scheme^ to generate smooth pseu- 
dopotentials, and remove the unphysical extremely non- 
local behaviour which results from the exact exchange in 
a controlled manner *± 

The resulting (tabulated) pseudopotentials are then 
represented in a convenient analytic form suitable for 
use in standard quantum chemistry codes, using a fit- 
ting procedure based on the method first described by 
Barthelat et al. 29 The version of this scheme employed 
here increases the efficiency and accuracy of the fitting 
procedure, and preserves the smoothness of the Troullier- 
Martins pseudopotentials. 

An analysis of the performance of our parameterized 
pseudopotentials for a number of test cases reveals that 
they reproduce the original AE atomic results well in 
most cases, and perform well in comparison with pseu- 
dopotentials published by previous authors which possess 
singularities at the nuclei. 

The pseudopotentials we have generated are finite at 
the nucleus and are smoother than other HF based pseu- 
dopotentials, which should aid the convergence of meth- 
ods employing Gaussian basis sets. These pseudopoten- 
tials are appropriate for use within non-relativistic theo- 
ries of the valence electrons, such as HF theory or corre- 
lated wave function methods, including DMC. 



Acknowledgments 

We would like to thank Dr. Y. Lee, whose Ph.Di^l 
work provided the algorithm which has evolved into the 
work presented here. Financial support was provided by 
the Engineering and Physical Sciences Research Council 
(EPSRC), UK. 



* Electronic address: jrt32@cam.ac.uk 

1 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 
(1980). 



2 W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Rajagopal, 
Rev. Mod. Phys. 73, 33 (2001). 

3 D. M. Ceperley, J. Stat. Phys. 43, 815 (1986). 



9 



B. L. Hammond, W. A. Lester, Jr., and P. J. Reynolds, 
Monte Carlo Methods in Ab Initio Quantum Chemistry 
(World Scientific, Singapore, 1994). 

C. W. Greeff and W. A. Lester, Jr., J. Chem. Phys. 109, 
1607 (1998). 

W. Miiller, J. Flesch, and W. Meyer, J. Chem. Phys. 80, 
3297 (1984). 

E. L. Shirley and R. M. Martin, Phys. Rev. B 47, 15413 
(1993). 

Y. Lee and R. J. Needs, Phys. Rev. B 67, 035121 (2003). 
L. Kleinman, Phys. Rev. B 21, 2630 (1980). 
G. B. Bachelet, D. R. Hamann, and M. Schliiter, Phys. 
Rev. B 26, 4199 (1982). 

P. J. Hay and W. R. Wadt, J. Chem. Phys. 82, 284 (1985); 
ibid. 82, 270 (1985); ibid. 82, 299 (1985). 
P. A. Christiansen, Y. S. Lee, and K. S. Pitzer, J. Chem. 
Phys. 71, 4445 (1979). 

L. A. Fefee, S. A. Wildman, G. A. DiLabio, T. M. Mof- 
fett, Jr., J. C. Peploski, and P. A. Christiansen, Clark- 
son University Relativistic Effective Potential Database, 
|http : // www, clarkson . edu/ ~pac/reps . html] 
M. Dolg, U. Wedig, H. Stoll, and H. Preuss, J. Chem. 
Phys. 86, 866 (1987); J. Sheu, S. Lee, and M. Dolg, J. 
Chin. Chem. Soc. 50, 583 (2003). 

M. Dolg and H. Stoll, Pseudopoten- 
tials of the Stuttgart/ Cologne group, 



V. R. Saunders, R. Dovesi, C. Roetti, M. Causa, N. M. 
Harrison, R. Orlando, and C. M. Zicovich- Wilson, CRYS- 
TAL98 User's Manual, University of Torino, Torino (1998). 

8 Gaussian 03, Revision C.02, M. J. Frisch, G. W. Trucks, 
H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheese- 
man, J. A. Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. 
Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, 

B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Pe- 
tersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. 
Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, 
O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. 
Hratchian, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, 
R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, 
R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. 
Morokuma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. 
G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, 
O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, 
J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clif- 
ford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, 
P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, 
M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challa- 
combe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, 

C. Gonzalez, and J. A. Pople, Gaussian, Inc., Wallingford 
CT, 2004. 

9 J. C. Barthelat, Ph. Durand, and A. Serafini, Mol. Phys. 
33. 159 (1977). 

http://www.theochem.uni-stuttgart.de/pseudopotentiale/inde^!fciGhtfigd]-Mann, H. Stoll, and H. Preuss, Mol. Phys. 65, 1321 

(1988). 



I. Ovcharenko, A. Aspary-Guzik, and W. A. Lester, Jr., J. 
Chem. Phys. 114, 7790 (2001). 

J. R. Trail and R. J. Needs, J. Chem. Phys., 122, 014112 
(2005). 

Y. S. Lee, W. C. Ermler, and K. S. Pitzer, J. Chem. Phys. 
67, 5861 (1977). 

C. Froese Fischer, T. Brage, and P. Jonsson, Computa- 
tional Atomic Structure: An MCHF approach (Institute of 
Physics Publishing, Bristol, 1997). 

G. B. Bachelet and M. Schliiter, Phys. Rev. B 25, 2103 
(1982). 

C. Frose Fischer, Comput. Phys. Commun. 43, 355 (1987). 
N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

A. L. Ankudinov, S. I. Zabinsky, and J. J. Rehr, Comput. 
Phys. Commun. 98, 359 (1996). 

C. Frose Fischer, Comput. Phys. Commun. 98, 255 (1996). 
L. Visscher and K. G. Dyall, At. Data Nucl. Data Tables 
67, 207 (1997). 

P. R. Librelon and F. E. Jorge, Int. J. Quantum Chem. 
95, 190 (2003). 



M. Dolg, H. Stoll, and H. Preuss, J. Chem. Phys. 90, 1730 
(1989); W. Kiichle, M. Dolg, H. Stoll, and H. Preuss, ibid. 
100, 7535 (1994). 

L. F. Pacios and P. A. Christiansen, J. Chem. Phys. 82, 
2664 (1985). 

M. M. Hurley, L. F. Pacios, P. A. Christiansen, R. B. Ross, 
and W. C. Ermler, J. Chem. Phys. 84, 6840 (1986). 
R. B. Ross, W. C. Ermler, and P. A. Christiansen, J. Chem. 
Phys. 84, 3297 (1986). 

J. P. Desclaux, Comput. Phys. Commun. 9, 31 (1975). 
See EPAPS Document No. E-JCPSA6-122-304518 for 
the parameters of our pseudopotentials and SO pseu- 
dopotentials. A direct link to this document may be 
found in the online article's HTML reference section. The 
document may also be reached via the EPAPS home- 
page (http://www.aip.org/pubservs/epaps.html) or from 
ftp.aip.org in the directory /epaps. See the EPAPS home- 
page for more information. 

Y. Lee, Ph.D. thesis, University of Cambridge, 2002. 



10 











Aei 


Ae-2 


At3 


Af4 


/\ffot 


Atom 


I 










xl0~ 


■5 




Li 


s 


-0.19634 


-0.19634 


0.0 


0.0 


-0.0 


0.0 


0.0 


Be 


s 


-0.30932 


-0.30933 


-0.1 


0.0 


-0.0 


-0.4 


-0.4 


n 


s 


-0.71212 


n 71 91 3 

u. f ±z ±0 


Q 
u. y 


n n 

u.u 


n n 

u.u 


n 4 


u.u 




P 


-0.40867 


-0.40869 


0.3 


0.1 


-2.2 


-0.5 


-2.3 


Si 


s 


u . ij^-'j ^ ^ 




n 9 


n n 

u.u 


n 1 

U. -L 


u.o 


n 9 

u.^ 






-0.27967 


-0.27968 


0.0 


0.3 


-1.1 


-0.6 


-1.3 


Ti 


s 


f) 99^{fiQ 
-u.^^ouy 


f) 99^7 


U.U 


-fl 


1 ^ 
-L .U 


1 n 7 


1 1 7 




d 


-0.39560 


-0.39542 


5.3 


0.1 


6.0 


6.1 


17.6 


(~V 


s 


-0.21104 


-0.21127 


-1.4 


n 1 

U. -L 


1 3 ^ 


7 Q 
- i . y 


-22.7 




d 


-0.31518 


-0.31494 


2.8 


0.6 


17.4 


2.8 


23.6 




s 




f) 98^7^ 

U.^Ut) ( 


1 ^ 
± .u 


n 1 

U. -L 


1.1 


n 7 

U. I 


n 4 




d 


-0.59311 


-0.59295 


8.7 


1.1 


3.3 


3.1 


16.3 


Br 




-1 09099 


-1.02101 


-0.8 


-1.5 


0.0 


0.3 


-1.9 




P 


-0.45616 


-0.45600 


-0.1 


16.2 


0.0 


-0.1 


16.0 


Mo 


s 


-0.21217 


fl 91 970 


U.U 


1.2 


-41.2 


1 1 Q 
- ± ± . y 


-52.4 




d 


-0.29622 


-0.29585 


0.5 


5.5 


36.8 


-5.5 


37.3 


As? 




-0.23716 


-0.23712 


-1.0 


3.9 


0.0 


1.6 


4.5 




d 


-0.51128 


-0.51108 


0.8 


21.1 


0.0 


-2.8 


19.2 


Sb 


s 


-0.62683 


-0.62689 


-0.3 


-1.2 


-5.3 


0.5 


-6.4 




P 


-0.30062 


-0.29996 


0.0 


54.7 


13.4 


-2.3 


65.8 






s- block 


0.1 


0.0 


0.0 


0.2 


0.3 


Mean Ae 




p- block 


0.8 


9.4 


1.2 


0.6 


11.1 






d- block 


1.9 


19.7 


7.4 


7.3 


29.0 



TABLE I: Eigenvalues (a.u.) for the neutral ground states of 
selected atoms. Eigenvalues obtained from HF pseudopoten- 
tial calculations with jj coupling are compared with AE DF 
eigenvalues. All eigenvalues shown are j-weighted averages. 
Errors are due to: 

Aei - transforming the Dirac equation to a Schrodinger equa- 
tion in the valence region. 

Ae2 - higher order effects due to excluding SO coupling. 
Ae3 - localization. 
Ae4 - parameterisation. 

Aei 0t = £; J — ef E = Yli A £ i is the total error. 
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TABLE II: Eigenvalues (a.u.) for the neutral ground states 
of selected atoms. The eigenvalues obtained from AE DF cal- 
culations are compared with those from HF pseudopotential 
calculations with jj coupling using our soft parameterized 
pseudopotentials, and with those from the singular parame- 
terized pseudopotentials of Pacios and Christianse n 12 ' 13 (PC). 
All eigenvalues shown are j- weighted averages. 
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AE HF 

Atom Configuration with RC Our PC 

Li 2~FfS] 0.00000 0.00000 0.00000 

2p 1 [ 2 P] 0.06767 0.06769 0.06782 

M 1 [ 2 D] 0.14076 0.14078 0.14089 

Be 2s 2 1 1 S] 0.00000 0.00000 0.00000 

2s 1 2p 1 [ 3 P] 0.06157 0.05951 0.06076 

2s 1 3d 1 [ 3 D] 0.23892 0.23910 0.23942 

C 2s 2 2p 2 [ 3 P] 0.00000 0.00000 0.00000 

2s 1 2p 3 [ 5 S] 0.08982 0.08461 0.08663 

2s 2 2p 1 3d 1 [ 3 P] 0.34003 0.34060 0.34102 

Si 3s 2 3p 2 [ 3 P] 0.00000 0.00000 0.00000 

3s 1 3p 3 [ 5 S\ 0.09313 0.09076 0.09243 

3s 2 3p 1 3d 1 [ 3 F] 0.21425 0.21458 0.21390 

Ti 3d 2 4s 2 [ 3 F] 0.00000 0.00000 0.00000 

3d 3 4s 1 [ s F] 0.02505 0.02945 0.02957 

3d 4 [ 5 D] 0.16459 0.18113 0.18015 

Cr 3d 5 4s 1 ! 7 S] 0.00000 0.00000 0.00000 

3d 4 4s 2 [ 5 D] 0.03914 0.04882 0.04820 

3d 6 [ 5 D] 0.26201 0.27226 0.26986 

Fe 3d 6 4s 2 [ 5 D] 0.00000 0.00000 0.00000 

3d 7 4s 1 [ 6 F] 0.07555 0.07984 0.08006 

3d 8 [ 3 F] 0.28933 0.30548 0.30477 

Br 4s 2 4p 5 [ 2 P] 0.00000 0.00000 0.00000 

4s 1 4p 6 [ 2 S'] 0.55464 0.55820 0.55803 

4s 2 4p 4 4d 1 [ 4 F] 0.33407 0.33376 0.33440 

Mo 4d 5 5s 1 [ 7 S] 0.00000 0.00000 0.00000 

4d 4 5s 2 [ 5 L>] 0.08763 0.10114 0.09298 

4d 6 [ 5 D] 0.15962 0.17646 0.18216 

Ag 4d 10 5s 1 [ 2 S] 0.00000 0.00000 0.00000 

4d 9 5s 2 [ 2 D] 0.14411 0.14888 0.15060 

4d 10 5p 1 [ 2 P] 0.11066 0.10873 0.11195 

Sb 5s 2 5p 3 [ 4 S'] 0.00000 0.00000 0.00000 

5s 1 5p 4 [ 4 P] 0.27820 0.28326 0.28370 

5p 5 [ 2 P] 0.76145 0.78007 0.78147 



Mean error- 


0.00513 


0.00539 


Mean absolute error 


0.00621 


0.00585 



TABLE III: Excitation energies (a.u.) for selected atoms. 
Excitation energies obtained from AE HF calculations with 
relativistic corrections (RC), are compared with those from 
HF pseudopotential calculations with LS coupling. Results 
are shown for both our soft parameterized pseudopotentials, 
and for the singular parameterized pseudopotentials of Pacios 
and Christiansen 1 ^ 3 - (PC). 
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0.19634 


0.19646 


Be 


0.29566 


0.29584 


0.29618 


C 


0.39625 


0.40032 


0.39724 


Si 


0.28083 


0.28743 


0.28083 


Ti 


0.20374 


0.20044 


0.20115 


Cr 


0.22054 


0.22849 


0.22757 


Fe 


0.23288 


0.22623 


0.22792 


Br 


0.39481 


0.39473 


0.39530 


Mo 


0.22801 


0.24066 


0.24459 


Ag 


0.23088 


0.23353 


0.23274 


Sb 


0.31787 


0.31923 


0.31852 


Mean 


error 


0.00231 


0.00188 


Mean absolute error 


0.00414 


0.00325 



TABLE IV: Comparison of first ionization potentials obtained 
from AE calculations {E io E ), from HF calculations using our 
soft parameterized pseudopotentials (Ei 0n ), and for the sin- 
gular parameterized pseudopotentials of Pacios and Chris- 
tianse n 12 ' 13 (Ef ^) (LS coupling is used throughout). 



Atom Configuration 


AE DF 


Our 


PC 


Li 


P l 


0.6 


0.8 


0.9 


Be+ 


P 1 


8.9 


8.5 


9.3 


B 


P 1 


20.6 


20.0 


21.8 


C+ 


P 1 


79.3 


75.8 


84.5 


N 2 + 


P 1 


209.8 


197.9 


219.9 


G 3+ 


P 1 


452.3 


454.4 


498.9 


F 


P' 


410.5 


402.4 


438.9 


Ne+ 


P' 


815.4 


820.5 


892.0 


Al 


P 1 


120.0 


116.1 


114.8 


Si+ 


P 1 


308.7 


293.9 


292.2 


p2+ 


P 1 


599.3 


556.1 


553.9 


S 3+ 


P 1 


1019.4 


1020.0 


1013.5 


CI 


P' 


920.6 


903.9 


910.3 


Ar+ 


P' 


1494.1 


1472.6 


1487.4 


Sc 


d 1 


165.7 


169.2 


191.0 


Ti+ 


d 1 


370.1 


410.7 


458.4 


Ti 3 + 


d 1 


405.4 


440.6 


487.3 


Cu 


d 9 


2122.9 


2110.5 


2292.9 


Zn+ 


rf" 


2820.0 


2814.1 


3036.3 


Ga 


P 1 


811.8 


756.6 


772.0 


Ge+ 


P 1 


1764.0 


1681.9 


1710.4 


Br 


P' 


3744.5 


3777.9 


3867.3 


Kr+ 


P 5 


5458.8 


5944.1 


5431.0 


Mean % error 

Mean absolute % error 


0.4 
4.8 


5.7 
8.2 



TABLE V: Spin-orbit splitting (cm - ) for selected atoms. 
Splitting obtained from all-electron DF calculations (with no 
relativistic correction) are compared with those resulting from 
jj-coupled HF calculations using parameterized AREP and 
SO pseudopotentials. Results are given for both our soft pa- 
rameterized AREP and SO pseudopotentials, and those of 
Christiansen et o/.— i 3 -^ (PC). 



